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1 Introduction 

An in-shutfle of a deck of cards is done by cutting the deck into two equal halves 
and interleaving them perfectly, with the first card of the second half being the 
first card of the shuffled deck. In case of 2n cards, this kind of a shuffle gives 
rise to a permutation (called the in-shuffle permutation of order 2n) given by: 

i -> 2i mod (2n + l),i G {1, 2, . . . , 2n} 

In [EM00 it has been shown that the problem of merging lists in-place can 
be reduced to performing in-shuffles in-place. In the same paper, the authors 
provide a linear time, in-place algorithm for performing the in-shuffle, when the 
list is represented as an array, indexed 1 to 2n. In this paper we provide a 
simpler algorithm which permits an easy implementation. 

2 Preliminaries 

Any permutation is a set of disjoint cycles. For e.g. the in-shuffle of 6 elements 
is composed of two cycles: (1, 2, 4) and (3, 6, 5). It is easy to move the elements 
of a single cycle using just an extra location, by a "cycle leader" algorithm, 
[FMP95J. This algorithm proceeds by repeatedly making a space in the list, 
computing the index of element that belongs to that space and moving that 
element, so creating a new space. 

For any permutation, we can proceed by applying the cycle leader algorithm 
to each cycle in turn in order to realise the whole permutation. If, for a permu- 
tation, we can compute the location of an element of the 'next' cycle in constant 
time and minimal space, we could realise the whole permutation using a linear 
time and an in-place algorithm. 

We show that, when 2n is of the form 3 fe — 1, we can easily determine the 
cycles of the in-shuffle permutation of order 2n. We will need the following 
theorem from number theory: 

Theorem 1 If p is an odd prime and g is a primitive root of p 2 , then g is a 
primitive root of p k for any k> 1. 



A proof of this theorem can be found in [NarOOl p 20-21]. 

It can be easily seen that 2 is a primitive root of 9. From the above theorem 
it follows that 2 is also a primitive root of 3 fc for any k> 1. This implies that 
the group (Z/i k )* is cyclic with 2 being its generator. 

Now let us analyse the cycles of an in-shuffle permutation when 2n = 3 —1. 

The cycle containing 1 is nothing but the group (Z/3 k )*, which consists of 
all numbers relatively prime to 3 k and less than it. 

Let 1 < s < k. Consider the cycle containing 3 5 . Every number in this cycle 
is of the form 3 S 2* (modulo 3 fc ) for 1 < t < <fi(3 k ) (where ip is the Euler-totient 
function). Since 2 is a generator of (Z/3 k )*, this cycle contains exactly the 
numbers less than 3 fe which are divisible by 3 s but not by any higher power of 
3. 

This means that in an in-shuffle permutation of order 3 fe — 1, we have exactly 
k cycles with 1, 3, 3 2 , . . . , 3 fe_1 each belonging to a different cycle. Thus for these 
permutations, it becomes easy to pick the 'next' cycle in order to apply the cycle 
leader algorithm. Note that the length of the cycle containing 3 s is (p(3 k )/3 s , 
which helps us implement the cycle leader algorithm more efficiently. 

We present the algorithm for the general case of an in-shuffle permutation 
of order 2n in the next section. 

3 Algorithm 

We start with an observation. Suppose for each 2n we can easily find a 2m (= 
£l(n) and < 2n) for which we have a linear time, in-place, in-shuffle algorithm, 
we then have a linear time, in-place, in-shuffle algorithm for 2n. This is because 
we can shuffle the array 

3 ®2 3 • • • 3 ) ) • • • 3 ) 3 • • • 3 ^n+m 3 3 • • • 3 ^2n] 

to look like 

ai, Q.2, . • • , a m , a„+i, . . . , a„ +m , a m+ i . . . , a„, a„+ m +i . . . , d2 n ] 

in linear time and in-place by doing a right cyclic shift of the elements a m +i, . . . , a„ + . 
by a distance m. We can achieve this by reversing the sub-array [m+1, . . . , n+m] 
followed by a reversal of the sub- arrays [m + 1 , . . . , n] and [n + 1 , . . . , n + m] . We 
can now perform an in-shuffle of the first 2m elements and recursively perform 
an in-shuffle of the remaining elements. 

Thus we have the following In-shuffle Algorithm: 

Input: An array A[l,..., 2n] 

Step 1. Find a 2m = 3 fc - 1 such that 3 fc < In < 3 fc+1 

Step 2. Do a right cyclic shift of A[m + 1, . . . , n + m] by a distance m 

Step 3. For each i G {0, 1, . . . , k — 1}, starting at 3*, do the cycle leader algo- 
rithm for the in-shuffle permutation of order 2m 

Step 4. Recursively do the in-shuffle algorithm on A[2m + 1, . . . , 2n]. 
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4 Space and Time Complexity 



Let T(n) be the time taken on an input array of size n. Let us look at the time 
and space bounds for each of the steps. 

Step 1 can be done in-place in O(logn) time. 

Step 2 can be done in-place in 0(n) time. 

Step 3 can be done in-place in O(m) time. 

Step 4 is a tail recursive call and hence can be implemented in constant space. 
The time taken by this call is T(2(n — to)). 

Since to = fi(n), it follows that the total time taken by this algorithm is 0(n). 

Thus we have a linear time, in-place algorithm for the in-shufne of a list of 
2n elements. 

5 Conclusion 

We have described a linear time, in-place, in-shufne algorithm which I believe 
is interesting, simple and easy to implement. It is not obvious whether the 
approach used (using a p such that k is a primitive root of p 2 ) in this paper can 
be generalised for fc-way shuffles, but it is probable that we can use this approach 
for fixed squarefree k, [HEA86] . For small k, this approach can definitely be 
used. 
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